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I . INTRODUCTION 



A. ARMT OPERATIONAL CONCEPT 

Conventional artillery uses unguided projectiles to 
destroy a variety of targets including armored vehicles. 
Ballistic information and professional experience are used 
to direct fire onto a target with little chance of a direct 
hit in the first few rounds. A forward observer provides 
feedback to the guns for correction until the target is 
destroyed with subsequent rounds . General ly this approach is 
limited to the 30 kilometer range and often point target 
locations are not known. 

The ARMT projectile, formerly called AIFS, changes the 
conventional concept by providing a fire and forget , boosted 
and terminally guided projectile with a highly predictable 
trajectory . Conventional heat seekers and target identifica- 
tion systems are short range and require a blunt nose which 
has high drag. Active radar is not practical . Long ranges of 
60 kilometers are achieved by using solid fuel ramjet 
propulsion with an axisymmetric air inlet spike. An 
accurately predicted trajectory is achieved by ensuring 
thrust equals drag to produce zero axial acceleration and 
create a pseudo- vacuum trajectory during powered flight. At 
a preset time (near ideal condition ramjet burnout) the 
inlet spike is discarded and canards are deployed for 
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terminal guidance control. A passive or reflected target 
energy system will be used for terminal guidance. Research 
is ongoing and sponsored by DARPA. Projectile 
characteristics are at Appendix A. 




Figure 1-1 ARMT Trajectory Concept 
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B. THRUST DRAG BALANCE CONCEPT 



Balanced thrust drag implies that the thrust equals the 
drag vectors during powered flight. The aim of the balance 
is to negate meteorological variations throughout the 
trajectory by reducing all force vectors to zero except 
gravity. The resulting pseudo- vacuum trajectory is highly 
predictable. The balance is accomplished by measuring axial 
acceleration and maintaining it at zero. Appendix B explores 
this concept in detail. The potential errors include 
transients due to wind induced reactions and measurement 
errors in the technical application of the balance in 
flight. Further, any asymmetries such as aerodynamic 
misalignment' or thrust will produce a lateral force which 
must be offset by a small roll rate(sixth degree of 
freedom) . This allows longitudinal forces to be the 
significant forces . 

The concept of thrust drag balance is best demonstrated 
in Figures 1-2 and 1-3. In Figure 1-2 the centre of gravity 
of the projectile is flying at an angle of attack (alpha) 
with a velocity vector (Vcg) relative to a fixed earth 
reference frame. The aerodynamic forces are lift(L) and 
drag(D). Assuming static static stability and neglecting 
transient dynamics, the angle of attack becomes zero, which 
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reduces lift to zero. If thrust equals drag then only the 
gravitational force exists. 

If a crosswind ( Vw ) is applied, the projectile cocks 
into the wind because of the aerodynamic force acting on the 
centre of pressure. Adjusting the thrust and again assuming 
static stability the drag can be offset as shown in Figure 
1-3 by applying the correct thrust. For future reference the 
component of drag normal to the plane containing Vcg is Dw . 




Figure 1-2 Projectile Dynamics Thrust Equals Drag 



Figure 1-3 Projectile Dynamics With Crosswind 



If a tailwind or headwind is applied, thrust is adjusted 
to compensate in the same manner. 

The dynamics or transient effects due to wind which 
create damped oscillatory motion have a cumulative net lift 
effect< as shown in Figure 1-4) as angle of attack is 
proportional to lift. The area under the curve in Figure 1-4 



does not average to zero due to the damping. 



For future 



reference this lift due to wind will be referred to as Lw . 




Figure 1-4 Angle Of Attack Versus Time 



Calculations integrating alpha(t) twice to demonstrate the 
above in the yaw or pitch plane are at Appendix B. 

C. THESIS SCOPE 

It is emphasized that the scope of this thesis is to 
analyze dispersion from the pseudo- vacuum trajectory at 
terminal guidance time for the eight inch ARMT projectile. In 
colloquial terms the "basket " from which the projectile 
will be guided to target was determined. Further, this 
basket is based on thrust errors and wind effects for a 
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launch angle of 45 degrees and muzzle velocity of MACH 2.2 



at sea level. Linearized aerodynamics were 
calculations. No attempt was made to establish 
for projectile impact. 



used in all 
a "footprint" 



20 



II. TRAJECTORY ANALYSIS METHODOLOGY 



A. VACUUM TRAJECTORY ANALYSIS 

A short and simple BASIC program was written to 
establish the vacuum trajectory and fuel burn time of the 
projectile in a two dimensional scenario where thrust 
exactly equaled drag. The Standard U.S. Atmosphere was used 
and the drag coefficient linearized [Reference 1]. The 
program listing and baseline data are at Appendix D. The key 
parameters and equations used are listed at Appendices A and 
D . The baseline pseudo- vacuum trajectory is shown in Figure 
2 - 1 . 

B . THREE-DEGREES-OF-FREEDOM MODEL 

The computer simulation was established using the 
Interactive Ordinary Differential Equation (IODE) library 
program at NPGS. A listing of the three degrees of freedom 
equations is in the IODE specifications at Appendix E. The 
simulation is based on wind tunnel results from Norden 
Systems [References 2 and 31 and assumes a horizontal plane 
at 30000 feet a lti tude . This height was chosen as significant 
flight time is at altitude and the time of flight in this 
plane was near actual flight time. Primarily the model was 
used to study crosswind effects but resonance with the 
projectile short period was simulated. Different thrust 
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scenarios were also used. Figure 2-2 depicts the model in 
the X-Y plane. 



burnout time=101 secs 
burnout range=175100 feet 
burnout height=10990 feet 

”o °. " 




Figure 2-1 ARMT Pseudo- vacuum Trajectory 




Figure 2-2 Three Degrees Freedom Model (Yaw Plane) 



The crosswinds were varied in intensity and duration to 
analyze crossrange dispersion. As well thrust was varied as 
a function of drag by scaling the drag. Cmq predictions were 
approximately -190 to -205 by Norden Systems [Reference 2] 
and an analysis was done to determine the effects of varying 
this damping coefficient. To further study the projectile 
characteristics, trajectories were simulated with no thrust 
or mass attenuation with varying wind profiles. Figure 2-3 , 



shows the baseline three-degrees-of -freedom trajectory o 
thrust equals drag with no wind. 




Figure 2-3 Baseline 



Three- Degrees -Of -Freedom Trajectory 



c. 



FI VE- DEGREES -OF -FREEDOM MODEL 



Again the IODE program was used to solve a five degrees 
of freedom system of equations. The roll plane was neglected 
for the following reasons: 

-the program could not accept sufficient functions to 
define six degrees of freedom. 

-the fins are expected to damp out initial spin due to 
gun rifling (900 RPM) in the first few seconds of flight. 
Sufficient spin to average projectile asymmetry errors to 
zero but not enough to significantly affect dynamics will 
remain. (ie. gyroscopic effects and Magnus effects) 

-the initial spin damping dynamics are beyond the scope 
of this thesis. 

The projectile was assumed to be axisymmetric allowing for 
ease of calculations in the yaw and pitch planes. Figures 
2-2 and 2-4 depict the yaw and pitch planes for the 
model, respectively . 

A listing of the equations is provided in the IODE 
specifications at Appendix F. Trajectory simulations were 
run with varying wind profiles and thrust error scenarios. 
Terminal guidance time (Tc) was determined by the baseline 
five degrees of freedom trajectory shown at Figure 2-1 as 
101 seconds. This trajectory has no wind and thrust exactly 



equaled drag until discard of the air inlet spike 



and 



canard deployment . A1 1 simulations assumed a 45 degree launch 



angle and Mach 2.2 muzzle velocity. Dispersion from the 
pseudo- vacuum trajectory was measured at terminal guidance 
time including the following variables: 



-X=range at 101 seconds 
-Y=crossrange at 101 seconds 
-Z=height at 101 seconds 
-Tb=burnout time 

-Vcg=speed at time 101 seconds 



Dispersion is defined as the difference between value for 
perturbed trajectories and the baseline trajectory. 




u Vw 




(body) 



L \ 

\ \ 



Z 



\ 




x (body) 



X (inertial) 



Vcg 



o 



Figure 2-4 Fi ve-Degrees-Of -Freedom Pitch Plane Model 



D. SIX-DEGREES-OF-FREEDOM MODEL 

i 

The Ballistics Research Laboratories (Aberdeen) six 
degrees of freedom model HTRAJ was obtained and adapted for 
use on the NPGS IBM 370 system. The adapted listing is 
included as Appendix G. Program compilation was achieved but 
output data was not satisfactory due primarily to input data 
syntax problems. Further, the program requires extensive 
modification to be able to vary thrust as a function of drag 
with the ease of the previous models already described. Time 
was not available to correct these problems. Detailed 
documentation on the program HTRAJ is available through BRL . 
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III. TRAJECTORY RESULTS 



A. THREE DEGREES OF FREEDOM 

1 • Crosawind And Thrust Error Sensitivity 

Figures 3-1 and 3-2 show the effects of constant 
thrust errors as a percentage of drag with no wind. Recall 
that the three-degrees-of freedom trajectory is at constant 
altitude. The figures therefore represent a thrust error 
sensitivity in terms of burn time and range variations due 
to drag or thrust. The trajectories all pass nearly through 
the 218000 foot range mark at t=100 seconds, but the burnout 
ranges are signif icantly different due to the varying burn 
times resulting from projectile thrust . Figure 3-2 
demonstrates the obvious conclusion that the burnout time 
decreases as thrust increases. 

Figure 3-3 combines the effects of constant 
crosswinds and thrust errors to demonstrate resultant cross 
range dispersion . The constant wind velocity lines are very 
nearly linear, and the pattern is essentially a shift about 
the zero-zero point. If thrust is less than drag, the 
dispersion is in the direction of the wind vector. If thrust 
is greater than drag then the overcompensation produces a 
dispersion opposite to the wind vector. 



£8 



RANGE (FEET) *10 

85.5 26.5 27.5 28.5 29.5 30.5 




Figure 3-1 Thrust Error Sensitivity (Range) 

Three Degrees Of Freedom 
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BURNOUT TIME (SECONDS) 




THRUST ERROR (%DRAG) 



Figure 3-2 Thrust Error Sensitivity (Burntime) 

Three Degrees Of Freedom 
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DISPERSION CROSSRANGE (FEET) 




Figure 3-3 Crosswind Sensitivity Three Degrees Of 
Freedom 
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Of interest was the dispersion due to the same wind 
profiles without thrust or mass loss. This represents an 
absolute dispersion with a given wind profile. These results 
are in Table 3.1. 



Table 3.1 No Thrust Crosswind Effects 
Three Degrees of Freedom 



Thrust=0 t 
vw (ft/sec) 

0 

60 

90 

120 



125 seconds 

dispersion (feet) 
crossrange 
0 

195.76 

293.80 

392.02 



The data clearly shows that crosswinds have a 
minimal effect in the thrusted scenarios . The dispersion in 
the thrust-equals- drag case is due to the transients (Lw) 



and is negligible as well. This is illustrated at Table 3.2. 



Table 3.2 Dispersion Due to Net Lift 

from Transients 



Thrust =drag 



vw (feet/sec) 
0 

60 

90 

120 



dispersion (feet) 
crossrange 

0 

.018 

.027 

.036 



Angle of attack damped to virtually zero during the 
trajectories (see Figure E-3) but the absolute maximum 
values were also recorded . Figure 3-4 demonstrates the 
relation between angle of attack (alpha) and the wind 
velocity . Thrust errors had virtually no effect on the angle 



of attack magnitude. 



ALPHA MAX (RADIANS) 







Figure 3-4 Wind Effects On Maximum Angle Of Attack 
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2 . Damping Coefficient Effects 

The magnitude of Cmq was determined to be 
approximately -200 in the flight speed ranges of the ARMT 
projectile [Reference 4] - To determine if this coefficient 
significantly affected dispersion Cmq was varied from 0 to 
1.9 x 10^ and the dispersion recorded. Figure 3-5 shows the 
relation between Cmq and dispersion (crossrange) with a 
constant wind of 60 feet/second and time of 125 seconds. 

o 




Figure 3-5 Damping Coefficient Effects 



Careful examination of Figure 3-5 shows that dispe- 
rsion increases slightly for realistic values of Cmq and is 
therefore proportional to Cmq . Due to the small sensitivity 
of dispersion to Cmq for all practical purposes one can 
conclude dispersion is independent to Cmq. (The increased 
dispersion proportional to Cmq can be explained as the 
difference in the areas under a highly damped sinusoid 
versus a non-damped sinusoid is considerable.) 

3. Sinusoidal Wind Analysis 

To test the effects of phase angles and wave- 
lengths on the projectile dynamics a sinusoidal wind was 
input to test the model. The wind (Vw) was of the form Vo 
sin<2 PI Vo t/lambda + phase) where Vo was set at 60 
feet/second . Figures 3-6 and 3-7 show the trajectories for 
the sine and cosine winds to burnout at thrust equal to 90% 
drag . (Lambda=range and lambda=.5 range)0f note is the fact 
that the sine wind produces a continuous dispersion while 
the cosine wind does not when thrust does not equal 
drag. This discovery led to an expanded study of different 
phase and wavelength sinusoidal wind inputs for various 
thrust scenarios. As well a simplified model was constructed 
as shown at Figure 3-8 and calculations confirmed the three 
degrees of freedom model results . Figures 3-9 and 3-10 depict 
the results of the calculations. 
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DISPERSION (FEET) 

3.2 4.0 4.8 5.6 6.4 7.2 8.0 




Figure 3-6 Sinusoidal Wind Effects (Sine) 
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o 

\n 




Figure 3-7 Sinusoidal Wind Effects (Cosine) 
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Y 




X 



Figure 3-8 Simplified Three Degrees Of Freedom Model 



The simplified model was used to analytically 
evaluate the crossrange dispersion (y) in the cases of sine 
and cosine winds to confirm the validity of the three- 
degrees -of -f reedom model . Calculations f ol low : 

General equations of motion in the y direction from 
Figure 3-8 are: 

M dVy/dt =Dy which approximately equals D(Vw-Vy)/Vxo 
if Vxo=Vo=constant . Then transforming the equation we get: 
dVy/dt + D/(MVo)Vy = D/(MVo) Vw(t) 
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or in time differentiated notation: 



y + B y = B Vw<t) if D/(MVo> =B 



If a cosine wind Vw(t)=Vwo cos (2 PI Vo t/A> is applied 
and the differential equation solved we have: 



y < t ) /A =K 



.-fit - 



cos ( 2 PI Vo t/A> + 



DAsin<2 PI Vo t/A) 
2 PI m Vo2 



let K= 2 1/2 m Vo2 Vwo 

DA Vo 



1 + f4 PI DA V 

[l/2 m Vo2 J 



Figure 3-9 graphically depicts the results as t increases. A 
zero asymptote is observed and this conforms with model 
results . 

Now if a sine wind Vw(t)= Vwo sin (2 Pi Vo t/A) is applied 
and the differential equation solved: 



y(t)/A =K 



DA 



(l-cos(2 PI Vo t/A> ) - 



4 PI 1/2 m Vo2 

sin ( 2 PI VO t/TO+ 4 PI(l/2 mVp2 ) ( j - e ~ Bt ) 

DA 



Figure 3-10 graphically depicts these results as 

increases . The asymptote is 

K (4 PI 1/2 m Vq2) (i- e -Bt) 

DA 

and this conforms to the model also. 
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Figure 3-9 Cosine Wind Analysis Results 
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Sine Wind Analysis Results 
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Figure 3-10 



The constant wind case was also derived: 



Vw(t) = Vw and the solved differential 

equation is: 

y(t)= Vw t - Vw m Vo/D ext- 
end as t increases y approaches Vw t. 

Specific calculations were done to confirm 
dispersion results at burnout. The calculations were based 
on : 



y=dy/dx X = 1/Vo dy/dt X 
Vwo=60 ft/sec 
Vo =2190 ft/sec 
>=304000 ft 
t= >/Vo 

D = l/2 j> Vo2 A Cd (=226 lbf) 
Thrust=.9D therefore D=22.6 lbf 
11=225 lbm 



Taking the time derivative of the cosine wind yCt) 
dispersion formula derived earlier we get 



dy/dt=K 



-Be-Bt t 2 PI Vo si n ( 2 PI Vo t/A> ♦ 

A 

DX . 2 PI Vo cos (2 PI Vo t/"\> 

4 PI 1/2 mVo^ * "X 
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at t= X/Vo the equation becomes: 



dy/dt = K(-Be~B^/Vo + b) 

and substituting the values into the approximat ionn 
for y; y=0 ft. 



Now , 


performing the same operations on the sine 


wind dispersion 


equation : 


dy/dt=K 


D"X _ 2 PI Vo sin ( 2 PI Vo t/"X) 

4 PI 1/2 mVo2 

-2 PI Vo cos ( 2 PI Vo t) + ( 4 PI 1/2 m Vo2). 
(Be-Bt) 

■ 



and at t = ?VVo 



dy/dt 


=K -2 PI Vo -4 PI 1/2 m Vo2.Be“Bt 

X D\ J 



and y=20 ft approximately. 

Finally, taking the constant wind dispersion 
formula and taking its derivative : 



dy/dt 


= Vw- Vw e'Bt 


and y= 10 ft 


approximately. 


These results 


are consistent with the three-degrees-of- 



freedom model considering the initial approximations. 

To study the short period resonance phenomenom the 
phase and wavelength wind variations were documented in 
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detail for a condition of thrust equal to 90% drag. Figures 
3-11 to 3-13 show the crossrange dispersion effects from 
wavelengths of 5000 ( lambda = .016 range) to 
303000 ( lambda=range ) feet and the phase angle variations 
from 0 to 2 PI. At approximately 9850 feet the dispersion 
character ist ics were no longer symmetric but remained under 
one foot. This is demonstrated at Figure 3-14. A crossrange 
dispersion contour was constructed at Figure 3-15 but this 
did not emphasize the resonance phenomenon sufficiently due 
to linearized aerodynamics and small dispersion 
magnitudes . The parameters that varied most significantly 
near resonance were angle of attack and burn t ime . Indeed , it 
is clear that the projectile is unstable in this region with 
a maximum angle of attack of 1.21 radians(70 degrees). Due to 
the increased angle of attack, drag increased and burn time 
decreased in the order of seven seconds (5%) compared to the 
baseline trajectory. This clearly affected the range at 
burnout. Due to large angle of attack the linearized 
aerodynamics are not val id . Pro j ect i le performance at wind 
resonance is worthy of additional analysis.lt was decided to 
demonstrate resonance by plotting the maximum angle of 
attack (alpha) versus wavelength with different thrust/ drag 
scenarios . The resonance results are at Figure 3-16. The 
projectile natural frequency is estimated at 0.22 Hz. 
Resonance was confirmed in the 0.22 Hz regime for thrust 
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errors under 10% including varying Vwo.This was necessary as 



the system is non-linear and Vwo could have affected 
resonance -It was noted that the curve shifted and increased 
or decreased in magnitude slightly in each case. This is due 
to the non-linearity of the system equations and the 
variations in Vo during flight. The absolute case of no 
thrust showed resonance at 10000 feet with a maximum angle 
of attack of 1.25 radians(72 degrees). 
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CROSSRANGE DISPERSION FEET 

Figure 3-11 Dispersion Versus Wind Phase Angle(Large Lambda) 



PHI RADIANS 

0.000 0.785 1.570 2.355 3.140 3.925 4.710 5.495 6.280 




CROSSRANGE DISPERSION FEET 



Figure 3-12 Dispersion Versus Phase Angle 
(Medium Lambda) 
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. PHI RADIANS 

0.000 0.785 1.570 2.355 3.140 3.925 4.710 5.495 6.280 




Figure 3-13 Dispersion Versus Phase Angle (Smal 1 Lambda 
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PHI RADIANS 

0.000 0.785 1.570 2.355 3.140 3.925 4.710 5.495 6.280 




CROSSRANGE DISPERSION FEET 



Figure 3-14 Dispersion Versus Phase Angle 

Resonance Lambda 
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PHI RADIANS 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 




Figure 3-15 Phase Angle Versus Wavelength 
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ALPHA MAX (RADIANS) 

0.25 0.50 0.75 1.00 1.25 




WAVELENGTH (FEET) 



Figure 3-16 Detrimental Resonance Between Projectile 
Short Period and Spatial Frequency of 
Wind Profile 
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It is clear that if the wind profile experienced by 
the projectile has harmonic content near the natural 
frequency then the round will be unstable and/or have large 
dispersion at burnout- This is best illustrated by viewing 
the wind as VwCt) with reference to the projectile. By a 
Fourier transform of the wind profile resonance with 
projectile natural frequency (fn) will be evident in the 
frequency domain. The wind oscillations need only be velocity 
variations with the appropriate harmonic content . Figure 3- 
17, which illustrates this concept, shows fn coinciding with 
a maximum of VwCf). One would anticipate stability problems 
for this case. 




A 

V 



FOURIER TRANSFORM (sample) 




Figure 3-17 Fourier Transform Of Wind Profile Concept 



Calculations for the long period resonance or 
phugoid show that the phugoid is irrelevant as flight time 
is too short: 

T = 4 . 44* Vo /g= 301 • 8 seconds [Reference 31 

The results of the sinusoidal winds were also 
compared with the aerodynamic jump phenomenom outlined in 
BRL report 1077 [Reference 4 ] . At first glance the 
sinusoidal wind results appear consistent with aerojump, but 
detailed calculations show otherwise. The rate of change of 
initial angle of attack is critical to aerojump and is 
easily incorporated into the aerojump equation. The results 
for the sine wind are a function of wavelength and of the 
same form as the simplified model results but the cosine 
wind and constant wind both yielded a small aerojump 
proportional to the initial wind. The magnitudes of the 
aerojump dispersions were more consistent with the net lift 
transient effects. Further, continuous perturbations in the 
form of wind after time zero are not in the aerojump 
equations and for all practical purposes the phenomenom is 
not observed in the model. These results should have been 
apparent before calculations were performed. Aerojump 
equations are not appropriate for non-zero wind 
conditions . Clear ly , when wind is included the aerojump 
equations cannot predict dispersion. 



B. FIVE DEGREES OF FREEDOM 



As previously stated dispersion was measured from the 
ideal pseudo- vacuum trajectory at a fixed time for 
initiation of terminal guidance (Tc) 101 seconds . Thrust 
error sensitivity in the zero wind condition is presented in 
Figures 3-18 and 3-19. In these two figures it is clear that 
thrust errors can have a significant impact on the location 
of * the projectile at the Tc. Figure 3-18 indicates that as 
thrust increases the dispersion moves downrange as 
expected . Figure 3-19 shows that as thrust increases the 
height of the projectile at Tc decreases. Due to added 
thrust the projectile is located at a point further along 
the parabolic trajectory; hence the height is less. 

As well Figure 3-20 shows these effects on available 
burn time of the ramjet. For thrust equals drag burnout at 
101 seconds exactly equals the Tc.With less thrust, terminal 
guidance is initiated before all the propellant is burned. 
Conversely, when thrust exceeds drag, burnout occurs prior 
to initiation of terminal guidance . Assuming this fixed time 
for initiation of terminal guidance, the projectile may 
experience up to seven seconds of unpowered flight ( 7 % 
error) before the guidance is deployed. The period of 
unpowered flight with associated high drag will cause major 
dispersions as shown in subsequent figures. 
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RANGE (FEET) 

173000 173500 174000 174500 175000 175500 176000 176500 




Figure 3-18 Thrust Error Effects (Range) 
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HEIGHT (FEET) 

1070(10750 10800 10850 10900 10950 11000 11050 11100 11150 11200 11250 




Figure 3-19 Thrust Error Effects (Height) 
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BURN TIME (SECONDS) 



z> 




THRUST ERROR (%DRAG) 



Figure 3-20 Thrust Error Effects (Burn Time) 
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Crosswind sensitivity results were very much the same 
as the three degree simulation . In fact the results appear 
even more linear with the same rotation about the zero zero 
point as in Figure 3-21. The gross features of Figure 3-3 and 
3-21 are identical but Figure 3-3 shows dispersion at 
burnout whereas Figure 3-21 shows dispersion at Tc. 

The crossrange dispersion results of a no thrust 
projectile with constant crosswind profiles and with thrust 
equal to drag (transients effects) are at Tables 3.3 and 3.4 
respectively . 



Table 3.3 No Thrust Wind Effects 

Five Degrees 



Thrust =0 
Vwxy (feet /sec) 

0 

30 

60 

90 

120 



time=101 seconds 

dispersion (feet) 
crossrange 

0 

92.51 

185.11 

272.90 

370.97 
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DISPERSION CROSSRANGE (FEET) 




Figure 3-21 Crosswind Effects Five Degrees 
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Table 3.4 Dispersion Due to Net Lift 

from Transients 



Thrust=drag time=101 secs 



Vwxy ( feet /sec ) 



dispersion (feet) 



0 



0 



30 



.0202 



60 



.0413 



90 



.0643 



120 



.0919 



Burn time was not significantly affected by crosswinds. 

Headwind and tailwind effects on range, height and 
burntime in different thrust scenarios are shown in Figures 
3-22 to 3-26. It is clear that dispersion from the pseudo- 
vacuum trajectory is minimized when thrust nearly equals 
drag. It is also apparent that the overthrusted trajectories 
involve greater dispersion as burnout occurs prior to 
terminal guidance and unpowered flight occurs for some time 
dependent on the scenario . This same error is observed in the 
case of headwinds and is most clear in the Figures 3-22B to 
3-25B about the thrust equal to drag axis. (Note that Figure 
3-22 to 3-25 inclusive have A,B and C sub-figures which are 
simply enlargements of key portions of the respective main 
figure ) . 
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RANGE (FEET) 

17380074000 174200 174400 174600 174800 175000 175200 - 175400 175600 175800 176000 176200 




Figure 3-22 



Headwind Sensitivity (Range) 
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Figure 3-22A Headwind Sensitivity (Range) 
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Figure 3-22B Headwind Sensitivity (Range) 
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Figure 3-22C Headwind Sensitivity (Range) 
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RANGE (FEET) 

17380074000 174200 174400 174600 174800 175000 175200 175400 175600 175800 176000 176200 




Figure 3-23 Tailwind Sensitivity (Range) 
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RANGE (FEET) 

174200 174300 174400 174500 174600 




Figure 3-23A Tailwind Sensitivity (Range) 
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Figure 3-23B Tailwind Sensitivity (Range) 
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Figure 3-23C Tailwind Sensitivity (Range) 
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APPENDIX D PSEUDO-VACUUM TRAJECTORY PROGRAM 



10 DIM X(300), Y(300), RO(300),D(300),MDOT(300),M(300), V(300) 
1£ DIM CDO(300) , MA (300) 

15 DIM T (300) 

£0 MMIN=198 
30 ISP=950 
40 DELT=. 5 
50 M ( 1 ) =££5 

&0 THETA=. 785398 1 G34000006# 

70 V0=£45£ 

75 LPRINT "TIME", "X", "Y", "MASS", "SPEED" 

80 FOR 1=1 TO 300 

90 X ( I ) =VO*COS ( THETA ) *T < I ) 

100 Y ( I ) = (VO*SIN (THETA) *T ( I ) ) - ( 1 6. 087*T ( I ) -'£) 

110 V(I) = ( (VO*COS (THETA) ) '"'£+ ( VO*S I N ( THETA ) - 
3£. 174*T(I) )'•'£) •■'.5 
1£0 TEMP =518. 688-. 003566£*Y(I) 

150 RO ( I ) = (£116. ££*<1-(.003566£*Y(I> /518. 688) ) --5. £56 ) / 
(1716. 5*TEMP) 

155 IF Y( I) >=36089! THEN RO ( I ) =. 00 1 £7 1 

157 IF Y(I)> =36089! THEN TEMP=389. 99 

158 MA ( I ) =V ( I ) /SQR ( 1. 4*1716. 3*TEMP) 

159 CDO ( I ) =. 50 1 5—. 1 0£5*MA ( I ) 

160 D ( I ) =. £76*. 5*R0 ( I ) *. 349 1*V ( I ) '£ 

161 D ( I ) =D ( I ) /. £76*CD0 ( I ) 

170 MDOT ( 1 + 1 ) =D ( I ) / ISP 

180 M ( 1+1 ) =M ( I ) -MDOT ( I ) *DELT 

190 T ( 1+1 ) =T ( I ) +DELT 

195 LPRINT T(I),X(I),Y(I),M(I),V(I) 

£00 IF M ( 1+1 ) <=MM I N GOTO ££0 
£10 NEXT I 

££0 LPRINT "RAMJET BURNOUT BETWEEN T=" ; T ( I ) ; " AND" ; T ( 1+1 ) 



9 £ 



burning . Burnout will occur on the vacuum trajectory but not 
at the same point, as shown in Figure C-2. 



z 




Figure C-2 Burnout Variations 
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In Figure C-l wind (Vw) has been added for the curve #4 
in the x direction . The new equations are : 

1*1^4 = C-D x + Dw) + T x 

• # 

JT124 = -D z + T z -mg 

t 

where Dx + Dw= 1/2 <^> U x 2 C<j A 

and UxCwind relative to the projectile) =Vx+Vw. With no 
thrust curve 4 results. 

Finally, if thrust is added to equal drag the 
following equations are valid: 

mx5 = 0 or T x = D x + Dw 
mz5 = -mg or T z = D z 

The accelerometer ( s ) read the acceleration caused by any 
differences in the thrust and drag forces (zero in this 
case) . 

It is clear that curves 1 and 5 are equal (and 3 if 
thrust=drag) but the burntime of a ramjet will not be equal 
in each case due to different drag profiles. 

t = (mass f uel )( Isp ) /thrust 

Therefore headwinds will require more thrust and 
burntime will be shorter . Conversely tailwinds allow longer 
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case Ux = Vx and Uz=Vz.This trajectory is represented by curve 



2. An accelometer in the x direction measures Dx/m and Dz/m 
in the z direction. 

If thrust is added, the trajectory is as shown at curve 
3< thrust greater than drag). The key equations are then: 
ra*X3 = -D x + T x 
mzg = -D z + T z -mg 

Accelerometers would measure (Tx-Dx)/m and <Ty-Dy)/m 
respectively . If thrust equals drag the trajectory would be 
curve 1 . 

z 




Figure C-l Pseudo- Vacuum Trajectory Concept Curves 

i 
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APPENDIX C-PSEUDO-VACUUM TRAJECTORY CONCEPT 



The basic concepts behind the pseudo- vacuum are not 
always readily apparent . This appendix analyzes the basic 
equations of motion in the following sequence of cases to 
illustrates the concept : 

vacuum trajectory 
atmosphere trajectory 
thrusted trajectory 
wind influenced trajectory 
thrust equal drag trajectory 

In Figure C-l (curve 1) the applicable equations of motion 
are : 

mx*i = 0 

mz*i = -mg 

An accelometer on the projectile in the z direction would 
measure no acceleration as it is in free fall already. This 
is similar to a satellite orbit or someone on an elevator in 
free fall. 

When atmosphere is added the key equations are: 

m#2 = -D x = -1/2 j> U x 2 A Cd 

•• _ 

ntZ2 = “D z - mg 

where U is relative wind and D is atmospheric drag. In this 
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When t=0, the equation reduces to: 

-2 S w +l(2$w) = 0 which is expected if the initial 
condition is zero. (ie. dispersion at zero equals zero at 
time zero ) . 

As t approaches infinity the equation value approaches 
infinite dispersion, but the oscillations of alpha have 
damped out- Note that the first derivative, proportional to 
velocity, tends to a constant value: 

( w /(<$2 + w 2 ))/m as t approaches infinity which is 
also to be expected. 
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APPENDIX B TRANSIENT ANGLE OF ATTACK CALCULATIONS 



The following equation of motion applies to the dynamics 
of angle of attack (alpha) oscillations due to the torque 
equation described in the three-degrees-of - freedom model. 



my = Lw which is proportional to alpha(t) 
where angle of attack is of the form : 



-St 

o<(t)=e sin(wt) where S is the damping constant 
and w is the angular velocity. Therefore: 



■ft- 

■a 



my = \ \ ©<(t) dt dT 

o 

e“St sin(wt) dt dT 



o 



■I 



e~<£T (-$sin(wT) - w cos(wT)) + w 
<5 2 + w2 



dT 



wt ( <S 2 +w2) -2Sw + 

e“^t (£2 s in ( wt ) + 2 S wcos ( wt ) -w^sin ( wt 

<TS2 + w 2)2 
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impulse = 950 "seconds'' 

Mass Property Equations: 

CG= < 225x39 . 4 ) - <52 . Oxf uel burned )/< 225-fuel burned) in 

Izz= . 35- < . 0019xf uel burned) in slug-ft 

Ixx= Iyy = 18 . 8- < . 0423xf uel burned) in slug-ft 

Aerodynamic Coefficient Equations: 

Cd = Cdo +.5 x Cl x sin( ) 

Cdo= -.1025 x MACH + .5015 <1.8< MACH < 2.2) 

Cm=Cl (Xcg-Xcp) /d 
C 1 = C 1. x «=< 

Key Wind Tunnel Results (per degree): 

MACH= 1 . 8 

Cdo= . 317 Cmq - - 203 

9 < 0 2 4 6 8 10 

Cl 0 .18 .38 .82 .92 1.3 

Xcp 45.4 45.6 48.2 49.7 50.2 50.0 

MACH=2 . 2 

Cdo= . 276 Cmq= - 190 

<*■ 0 2 4 6 8 10 

Cl 0 .19 .40 .65 .95 1.3 

Xcp 42.9 42.9 44.6 45.1 46.0 44.6 

Cmq=d(Cm)/d (qd/2V ) 

Initial Muzzle Spin Rate 30 PI radians /sec 
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APPENDIX A-PROJECTILE CHARACTERISTICS 



GUIDANCE AND CONTROL 



OBTURATOR 




Figure A-l ARMT Projectile Baseline Design 



The following characteristics are obtained from 
References 2 and 6. 
length=72 inches 
diameter=8 inches 
launch weight=225 Ibm 

launch CG=3S.4 inches from spike apex 

fuel weight=27 lbm (CG at 52.0 inches from apex) 
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except where harmonic content near resonance is 
contained in the wind profile. 

-headwinds create greater dispersion as ramjet 
burnout occurs before terminal guidance is 

initiated . Crosswinds and tailwinds generally will 
have minimal effects. 

Recommendations for further study related to this 
thesis are: 

-Dynamics at launch including the initial spin 
damping, Magnus effects and aerodynamic jump. 

-Dynamics at terminal guidance including air inlet 
spike ejection and canard deployment. 

-Complete adaptation of BRL HTRAJ on the NPGS IBM 
system to use non-linear aerodynamics. 
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IV. CONCLUSIONS 



Generally this thesis confirms the viability of the 
pseudo-vacuum trajectory from a thrust equals drag 
pro jecti le . Indeed winds up to 120 feet per second and thrust 
errors up to five percent created dispersions of less than 
600 feet over a range of 175090 feet for the ARMT 
projectile . 

Specifically the following conclusions are made: 

-the drag due to wind (Dw) is the aerodynamic 
force which creates the most dispersion. 

-the net lift due to transients (Lw) is a minor 
aerodynamic force in creating dispersion. 

-the dispersion is proportional to the damping 
coefficient C Cmq ), although for practical values of Cmq 
dispersion can be treated as insensitive to Cmq. 

-the projectile resonance or short period is at 
approximately 0.22 Hz and is unstable in the 
resonance regime due to variable wind profiles. The phugoid 
is irrelevant. 

-thrust errors are the most significant errors 
contributing to dispersion. 

-winds have a minor effect on dispersion during 
powered flight when thrust approximates drag. 
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Figure 3-28 Three Dimensional Trajectory With Altitude 

Dependent Wind Profile 
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Wind effects during near thrust equals drag portions of 
the trajectory were very small. The critical factor 
affecting the dispersion will be the setting of the terminal 
guidance clock in the pro jecti le . This presents several 
options for 'setting initiation of terminal guidance: 

-Fixed time irrespective of trajectory (ie. 
elevation of gun). Burn time would be shorter at less than 45 
degree launch and the basket would be more sensitive to wind 
and thrust errors. As well the trajectory would change and so 
would the predicted basket. 

-Fixed time for a given tr a j ectory . Thi s requires 
analysis of the sensitivity to elevation of the gun and a 
manual Tc setting method for gunners. 

-Ramjet burnout. This would ensure that Tc occurred 
at an ••optimum" point almost exactly on the pseudo- vacuum 
trajectory but the location along it at Tc would be 
difficult to predict. 

-At given pressure for altitude [Reference 5 ].This 
would establish accurate baskets for given gun elevations 
assuming the altitude chosen was before ramjet burnout. 

Several trajectories were simulated with a varying wind 
profile as a function of altitude in the Standard 
Atmosphere . These are shown at Figure 3-28 as compared to the 
pseudo- vacuum trajectory . 



Velocities at time 101 seconds were measured and are 
shown at Figure 3-27. The total variation of extreme cases is 
from 2285 feet/sec to 2315 feet/sec, or only 30 feet/sec. 




Figure 3-27 Velocity Sensitivity 
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Deviations from the ideal trajectory were correlated 
and tabulated in Table 3.5. Results were conservatively 
linearized and are presented in that manner. Wind errors are 
read in addition to the thrust errors at the left of the 
table. For example: 

If thrust error=-5x 

crosswind=120 ft/sec 

predicted dispersion is X=-540 ft 

Y = 19 ft 

2=130 ft ;from ideal basket. 

Table 3.5 Dispersion Analysis Results 





Thrust 


Errors 


Winds 

continuous 120 


ft /sec 




Errors 


X 


2 


Cross 

Y 


X 


Head 

2 


Tail 

X 


2 


0 * 


0 


0 


. 1 


30 


7 


. 1 


. 1 


1 % 


110 


25 


4 


30 


15 


10 


1 


2 % 


220 


50 


8 


40 


20 


20 


2 


5 % 


540 


130 


19 


50 


32 


50 


10 


10 % 


1100 


345 


38 


90 


50 


100 


30 
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Figure 3-26 Tail/Headwind Effects (Burntime) 
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Figure 3-25C Tailwind Sensitivity (Height) 
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Figure 3-25B Tailwind Sensitivity (Height) 
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Figure 3-25A Tailwind Sensitivity (Height) 
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Figure 3-25 Tailwind Sensitivity (Height) 
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Figure 3-24C Headwind Sensitivity (Height) 




Figure 3-24B Headwind Sensitivity (Height) 
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Figure 3-24A Headwind Sensitivity (Height) 
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Figure 3-24 Headwind Sensitivity (Height) 
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REPRODUCED AT GOVERNMENT EXPE 



IME 


X 


Y 


MASS 


SPEED 


0 


0 


0 
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a 5 


S66. 913 


362. 391 
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2440. 65 


1 


1733. 3 3 


1717. 74 


224. 637 
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APPENDIX F-FIVE CEGREES OF FREECOM RESULTS 
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